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We present the results of a theoretical study of a four-level atomic system in vee + ladder con¬ 
figuration using a density matrix analysis. The absorption and dispersion profiles are derived for a 
weak probe field and for varying strengths of the two strong control fields. For specificity, we choose 
energy levels of ®^Rb, and present results for both stationary atoms and moving atoms in room 
temperature vapor. An electromagnetically induced absorption (EIA) peak with negative disper¬ 
sion is observed at zero probe detuning when the control fields have equal strengths, which switches 
to electromagnetically induced transparency (FIT) with positive dispersion (due to splitting of the 
EIA peak) when the control fields are unequal. There is significant linewidth narrowing in thermal 
vapor. 
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I. INTRODUCTION 


The phenomenon of electromagnetically induced trans¬ 
parency (EIT), where the absorption of a weak probe 
laser on one transition is reduced using a strong con¬ 
trol laser on an auxiliary transition, is well-studied 
in three-level systems, in all the three con:^urations 
namely lambda (A), vee (V), and ladder (S) [ll-Q. But 
the related phenomenon of electromagnetically induced 
absorption (EIA)—enhanced absorption of the probe 
laser—is possible only when there are at least four lev¬ 
els, so that additional control lasers can be used. EIA has 
been studied—both theoretically and experimentally—in 
N-type (A -I- V) four-level systems M. The other four- 
level systems, namely Y-type (5 -|- 5), inverted Y-type 
(A -h S), and tripod-type (A -|- A), only show enhanced 
EIT. 


In this work, we present theoretical analysis of a new 
kind of four-level system formed by the combination of V 
and 5 three-level systems, which shows EIA resonances. 
This configuration opens up the possibility of the experi¬ 
mental study of EIA in an important class of atoms called 
Rydberg atoms, because such atoms have a highly ex¬ 
cited atomic state allowing the formation of a 5-type 
system. Rydberg atoms have thus been used to study 
EIT phenomena 0i, and the present analysis allows 
the extension to the observation of EIA in such atoms. 


For application to a real system, we choose energy lev¬ 
els of ®^Rb, and present results for both stationary atoms 
and moving atoms in room temperature vapor. An EIA 
peak with negative dispersion is observed at zero probe 
detuning when the control fields have equal strengths. It 
switches to an EIT peak with positive dispersion (due 
to splitting of the EIA peak) when the control fields are 
unequal. Thus, the atomic medium can be switched from 
EIA and super-luminal light propagation to EIT and sub¬ 
luminal light propagation. 



FIG. 1. (Color online) Four-level system formed by the com¬ 
bination of three-level vee and ladder systems. 

II. THEORETICAL CONSIDERATIONS 

The four-level system considered in the present work 
is shown in Fig. [TJ The ground state |1) is coupled with 
state |2) by a weak probe field. Two strong control fields 
are present—control 1 between levels |1) and |3), and 
control 2 between levels |2) and |4). The control fields 
have detunings Ad and Ac 2 , while the frequency of the 
probe field is scanned with variable detuning Ap. The 
intensities in the various beams are given in terms of the 
respective Rabi frequencies, which are flp, Od, and Oc 2 . 

The total Hamiltonian after carrying out the rotating 
wave approximation (RWA) is: 

H=^[^p |1) (2| + Hd |1) (3| + Hd |2) (4|] + h.c. 

+h [Ap |2) (2| + Ad |3) {3| -b (Ap -I- Ad) |4) (4|](1) 

where h.c. is the Hermitian conjugate of the preceding 
off-diagonal terms. As expected, in the absence of the 
control 1 (or the control 2) field, the four-level system 
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reduces to a three-level ladder (or vee) system. The time 
evolution of the system is given by the standard Liouville 
equation for the density matrix p. The coupled density 
matrix equations for this four-level system are given in 
the appendix. 

The observable in our work is the response of atoms 
to the weak probe field. The susceptibility of probe field 
is determined by the coherence between levels |1) and 
12)— Pi 2 —and is defined as 0 


; 1, 7 ^lMl2p 

X = kpi 2 where k = ——— (2) 

ftSQi up 

Here N is the atom number density in the medium and 
Pi 2 is the dipole matrix element for the probe transi¬ 
tion. The dispersion of the probe held is proportional to 

I 


Re{pi 2 }, and its absorption is proportional to Im{pi 2 }. 
The group velocity of the probe held is given by 
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where c is the velocity of light and ujp is the probe fre¬ 
quency. The above equation shows that Vg is inversely 
proportional to the slope of probe dispersion—positive 
dispersion (EIT in the imaginary part) gives rise to sub¬ 
luminal propagation, while negative dispersion (EIA in 
the imaginary part) gives rise to super-luminal propaga¬ 
tion. 

The various density-matrix elements are solved follow¬ 
ing the procedure given in Ref. 0- The steady-state 
solution for pi 2 is given by 
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Also, the population difference between states |1) and |3) 
is given by [nj 


-1 



As expected, the above equation yields the solutions for 
three-level S and V-systems by considering that only the 
corresponding control helds are non-zero [Tl, [l^ . 

For specihcity, we have chosen energy levels of the ®^Rb 
atom—level |1) is the F = 1 hyperhne level of the 5Si/2 
ground state; level |2) is the F = 2 hyperhne level of the 
5 P 3/2 state; level |3) is the F = 2 hyperhne level of the 
5 Pi/ 2 state; and level |4) is the F = 1 hyperhne level 
of the 5 D 5/2 state. The corresponding decay rates are: 
Pi = 0 (because it is a ground state), r 2 = 27r x 6.1 
MHz, Ps = 27r X 5.9 MHz, and r 4 = 27r x 0.68 MHz. 
Since the probe held is weak, it ensures that P 22 = 0, 
and hence P 44 is also 0. The population cycles between 
|1) and 13) (closed transition) according to Eq. ([5]). The 
energy levels considered here do not allow the formation 
of the well-known diamond conhguration dealt with in 
references [l^ and [T^, because the 5 D 5/2 state does not 


couple to the 5Pi/2 intermediate state. 

The above analysis is also applicable to Rydberg atoms 
with minor changes in parameters for level |4). For ex- 
ample, level |4) can be the 44 D 5/2 state considered in 
Ref. with a decay rate of r 4 27r x 0.3 MHz. 


III. RESULTS AND DISCUSSION 
A. Stationary atoms 

We hrst consider the results for stationary 
atoms. The imaginary and the real parts of the 
solution—corresponding to absorption and dispersion, 
respectively—as a function of probe detuning are shown 
in Fig. [21 The two control helds are taken to be on 
resonance, i.e. Ad = Ac 2 = 0. Parts (a) and (b) show 
the results when only one control held is present, i.e. 
corresponding to ladder and V-type three-level systems 
respectively. As expected, probe absorption shows an 
EIT dip near line center, because absorption splits 
into an Autler-Townes doublet at the locations of the 
dressed states created by the control held, i.e. when 
Ap = ±Hc/2. The corresponding dispersion prohle has 
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a positive slope near zero detuning, and hence can be 
used for sub-luminal light propagation. 

Now look at what happens when both control fields 
are applied. The results for fixed value of flci = 6 r 2 
and varying VLc 2 are shown in part (c)“(f) of the fig¬ 
ure. When the two control fields are equal, the medium 
switches from transmitting to absorbing, and has an EIA 
peak at line center. The corresponding dispersion profile 
has negative slope near zero detuning, which can there¬ 
fore be used for super-luminal light propagation. On 
either side of this equality condition, the central EIA 
peak splits into two, and shows a transparency dip at 
line center. The corresponding dispersion profile has pos¬ 
itive slope. Therefore, the medium can be switched from 
transparency (slow light) to absorption (fast light) by 
changing the relative strengths of the two control fields. 

B. Thermal averaging for moving atoms 

We now consider the results of Doppler averaging for 
an atomic vapor at room temperature, which is the typi¬ 
cal condition for Rb atoms used for an experimental real¬ 
ization. Since the probe and control beams are taken to 
be co-propagating or counter-propagating, the relevant 
distribution is the onc-dimensional MaxwelDBoltzmann 
velocity distribution. For an atom of mass M at a tem¬ 
perature T and moving with a velocity u, this is given 
by 

f{v)dv = J ^ exp dv 

' V 27rfcBT ^ V ‘^ksT J 

where kg is the Boltzmann constant. If the velocity is 
along the beam direction, then the detuning of the field 
will change by the Doppler shift of ±kv, where k is the 
photon wavevector and the sign depends on the relative 
direction of the beam with respect to the velocity. The 
thermal averaging over the MaxwelDBoltzmann distribu¬ 
tion is done by assuming that the probe field and control 
field 1 are co-propagating, and control field 2 is counter- 
propagating. 

The absorption profiles after thermal averaging, shown 
in Fig. [3l confirm our earlier observation of linewidth 
reduction in Ref. [la. The main results are the same 
as that seen for stationary atoms—an EIT dip for a 
three-level system when one of the control fields has zero 
strength; an EIA peak when the two control fields are 
equal; and the EIA peak splitting into two when the con¬ 
trol fields have unequal strengths. The linewidth reduc¬ 
tion is because moving atoms fill in the transparency win¬ 
dow. This can be seen in the figure, where the profile for 
atoms moving only with positive velocities is shown with 
a dashed line while the profile for atoms moving only 
with negative velocities is shown with a dotted line— 
these profiles are calculated after Doppler averaging for 
one sign of velocity only. 

The same linewidth reduction is also seen for the dis¬ 
persion profiles, shown in Fig. [D As in the case of sta- 



FIG. 2 . (Color online) Calculated probe absorption (Im{pi2}, 
solid), and dispersion (Re{pi2}, dashed) versus probe detun¬ 
ing for (a) Del = 0, Dc2 = 6r2; (b) Dd = 6r2, Dc2 = 0; 
(c) Del = 6r2, Dc 2 = 4r2; (d) Dei = 6r2, De2 = 6r2; (e) 
Del = 6 r 2 , De2 = 8 r 2 ; (f) Del = 6 r 2 , De2 = 12 r 2 . 
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FIG. 3 . (Color online) Effect of Doppler averaging on probe 
absorption at room temperature, (a) Del = 0, Dc2 = 6r2; 
(b) Del = 6r2, Dc 2 = 0; (c) Dd = 6r2, Dc2 = 4r2; (d) Dci = 
6r2, Dc 2 = 6r2; (e) Dci = 6r2, Dc 2 = 8r2; (f) Dci = 6r2, 
Dc2 = 12 r2. Solid curves are for all atoms, dashed curves are 
for atoms moving with positive velocities, and dotted curves 
are for atoms moving with negative velocities. 


tionary atoms, the slope near line center changes from 
positive for EIT to negative for El A, demonstrating the 
potential of switching the speed of light propagation. 


IV. CONCLUSION 


In conclusion, we have investigated a four-level atomic 
system in vee -I- ladder configuration, with a weak probe 
field and two strong control fields. For application to a 
real system, we consider relevant low-lying energy lev¬ 
els in ®^Rb. The coupled density-matrix equations are 
solved in steady state to calculate the absorption and 
dispersion profiles of the probe beam. The calculation is 
done for both stationary atoms and for atoms moving in 
room temperature vapor. In both cases, probe response 
shows an EIA window with negative dispersion at line 
center when the control fields have equal strengths. The 
induced absorption peak splits when the control fields are 
unequal, resulting in an EIT window with positive disper¬ 
sion. For moving atoms, Doppler averaging leads to sig¬ 
nificant linewidth narrowing for both absorption and dis¬ 
persion profiles. These results can therefore be used for 
applications such as switching between sub- and super¬ 
luminal light propagation inside an atomic medium. 
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APPENDIX 


In our analysis, the coupled density matrix equations 
in rotating frame are calculated by using the definition 
of density matrix and the expression for Hamiltonian. 
The time evolution of system—by incorporating the de¬ 
cay of atoms from each level and repopulation from ex- 
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FIG. 4 . (Color online) Effect of Doppler averaging on probe 
dispersion at room temperature, (a) fid = 0, Dc2 = 6r2; 
(b) flci = 6r2, 0,c2 = 0 ; (c) fici = 6r2, = 4r2; (d) flci = 

6r2, IIc2 = 6r2; (e) fid = 6r2, flc2 = 8r2; (f) Dd = 6r2, 
r2c2 = 12 r2. Solid curves are for all atoms, dashed curves are 
for atoms moving with positive velocities, and dotted curves 
are for atoms moving with negative velocities. 


cited levels—is given by 


Pll = r 2 P 22 + r 3 P 33 + - (n*pi 2 — npP 2 l) 

+ 0 (^ciPi3 - ^ciPii ); 


P22 — —^2P22+^APA4:+-^{P‘pP21—^*pPl2) 

+ 2 (^*2^24 — ^c2PA2) ; 

P33 = -rspss + 2 (^ciP3i — ^ciPis ); 

P44 = —r4P44 + - (flc2P42 “ ^*c2P2a) 1 

Pl2 = f-^ + *Ap j pi2 + (pil - P22) 


Pl3 = 


“2 (^ciP32 - ^^* 2 ^ 14 ); 

r3 \ i 

—+ *Acl j pi 3 + 2 ^'=! (^11 ~ ^33) 


Pl4 = 


--j^pP23] 

■p 

-^ + *(Ap + Ac 2) ) Pl4 


+ 2 (^c2Pl2 — ^pP2A — ^cIPZa) ; 


+ *(Acl — Ap) I P23 


P 23 = 


+ 2 {P-clP2l — ^pPl3 — ^c2PA3) ; 

/r 2 +r 4 •a') *rA/ ^ 

P 24 = I-^- 1 -iAc 2 I P 24 + 2 ^^c 2 (P 22 — P 44 ) 


P3A = 


--Vi*pPiA-, 

r3 + r4 


+ *(Ap + Ac 2 — Aci) ) P34 


+ 2 (^c2P32 - ^*ciP1a) ; 


The positive decay terms in pa {i = 1, 2, 3,4) equations 
are due to the fact that there is an increase in population 
of lower levels due to the decay of atoms from higher 
levels. 
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